arXiv:1506.01473v2 [gr-qc] 2 Feb 2016 


Comparing Hamiltonians of a spinning test particle for different tetrad fields 

Daniela Kunst/’B Tomas Ledvinka,^’0 Georgios Lukes-Gerakopoulos,^’^’0 and Jonathan Seyrich^’B 

^ ZARM, University of Bremen, Am Fallturm, 28359 Bremen, Germany 
^Institute of Theoretieal Physies, Faculty of Mathematics and Physics, 

Charles University in Prague, 18000 Prague, Czech Republic 
^Institute of Theoretical Physics, Faculty of Mathematics and Physics, Charles University in Prague, Czech Republic 
^ Theoretical Physics Institute, University of Jena, 077f 3 Jena, Germany 
^Mathematisches Institut, Universitdt Tubingen, Auf der Morgenstelle, 72076 Tubingen, Germany 

This work is concerned with suitable choices of tetrad fields and coordinate systems for the 
Hamiltonian formalism of a spinning particle derived in [E. Barausse, E. Racine, and A. Buonanno, 

Phys. Rev. D 80, 104025 (2009)]. After demonstrating that with the originally proposed tetrad 
field the components of the total angular momentum are not preserved in the Schwarzschild limit, 
we analyze other hitherto proposed tetrad choices. Then, we introduce and thoroughly test two new 
tetrad fields in the horizon penetrating Kerr-Schild coordinates. Moreover, we show that for the 
Schwarzschild spacetime background the linearized in spin Hamiltonian corresponds to an integrable 
system, while for the Kerr spacetime we find chaos which suggests a nonintegrable system. 
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1. INTRODUCTION 

The motion of a spinning particle in the spacetime 
background of a black hole, particularly the Kerr space- 
time, is of great astrophysical interest. Namely, it ap¬ 
proximates the motion of a stellar compact object (e.g., 
a black hole) around a supermassive black hole. Such 
binary systems are expected to lie at the center of galax¬ 
ies (see, e.g., [l| and references therein) and to be good 
candidates for sources of gravitational radiation Q . 

Even though the equations describing the motion of a 
spinning particle in a curved spacetime have been pro¬ 
vided several decades ago by Mathisson Q and Papa- 
petrou 0 , many issues of this motion are still open. The 
problem lies in the fact that the Mathisson-Papapetrou 
(MP) equations are not a closed system of first order 
differential equations. Hence, a spin supplementary con¬ 
dition (SSG) is needed in order to close them. Several 
such SSGs have been proposed (see, e.g., d, @ for a 
review), each of which introduces a different reference 
frame. Physically, the ambiguity in the choice of a SSG 
is related to the fact that a spinning body cannot be 
treated as a point particle but must have a finite size 
in order to be prevented from rotating at superluminal 
speed 0. In particular, each SSG corresponds to an ob¬ 
server who sees the reference worldline fixed by the SSG 
as the center of mass of the extended body. Previous 
studies have shown that the choice of the SSG depends 
on the question one wants to investigate [Ml. 

Within the MP equations the motion of spinning test 
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particles in the Schwarzschild or Kerr spacetime have 
been investigated in [iM3, and several papers have 
been devoted to the investigation of the appearing chaotic 
motion [l6l - [l9l| . Beyond the pole-dipole approximation, 
the quadrupole moment of the test particle has already 
been taken into account [lO, [Hj . 

The dynamics of spinning test particles has not only 
been worked out in Lagrangian formalisms (MP equa¬ 
tions) [l^, but in Hamiltonian formalisms [23l - l^ as well. 
Hamiltonian dynamics has a long tradition in astronomy 
and a large number of problems there (e.g., perturbative 
problems or chaotic motion) are typically studied from 
a Hamiltonian perspective [29|. In general relativity the 
Hamiltonian formalisms have been applied, for example, 
in the framework of the canonical Arnowitt-Deser-Misner 
(ADM) formalism 1)^ and in the effective-one-body ap¬ 
proach (EOB) |3ll433l| . which studies the dynamics of 
spinning bodies of compact objects using the Hamilto¬ 
nian description of a one-body problem [3^. Due to the 
significance of a Hamiltonian approach the Hamiltonian 
description of spinning particles is important despite the 
fact that it mostly neglects terms quadratic in spin. 

In our previous work Q , we have compared the Tulczy- 
jew (T) SSG [ 3 ^ with the Newton-Wigner (NW) SSG 
as supplements to the MP equations. In a second step, 
we compared the MP equations supplemented by the NW 
SSG to the corresponding Hamilton equations derived in 
based on the same NW SSG. In this work we focus 
on the latter, i.e., on a canonical Hamiltonian formalism 
which should be equivalent to the MP equations up to 
the linear order of the test particle spin. 

In contrast to the T SSG, the NW SSG, which is used 
within the framework of the Hamiltonian formalism, does 
not provide a unique choice of reference frame. It rather 
dehnes an entire class of observers, each characterized by 
a different tetrad field. Thus, the Hamiltonian formalism 
proposed in depends on the choice of a reference basis 
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given by such a tetrad field. Each choice of a tetrad field 
basically determines the form and the properties of the 
resulting Hamiltonian function. The fact that tetrads 
providing certain frames of reference are involved in a 
definition of the spin variable can also be seen as a con¬ 
sequence of the fact that in the Hamiltonian description 
the spin is a vector with prescribed canonical relations to 
coordinates and momenta. Still, one might conclude that 
the tetrad dependence of the Hamiltonian description of 
the spinning particle is against covariance principles of 
general relativity. Yet, when we numerically solve equa¬ 
tions of motion we have to use some coordinates anyway. 
The involvement of tetrads simply means we use different 
coordinates for external and inner degrees of freedom of 
the spinning particle. As, e.g., Boyer-Lindquist coordi¬ 
nates are comfortable for solving equations of motion in 
Kerr geometry there may well be some other tetrad fields 
more suitable for the definition of the Hamiltonian spin. 

We discuss the advantages and the drawbacks of 
Hamiltonian functions arising from tetrad fields already 
proposed in [ 2 ^ [s^ . Then, we introduce two new tetrad 
fields in Kerr-Schild coordinates which yield Hamiltonian 
functions with desirable properties using both analytical 
and numerical analysis. Namely in order to have a good 
choice of a tetrad field, the corresponding Hamiltonian 
should reflect the symmetries of the background space- 
time, i.e., preserve the integrals of motion, and avoid any 
coordinate effects evoked by coordinate dependent tetrad 
basis vectors. For the above discussion, we focus on 
the Schwarzschild limit and show that the well behaving 
Hamiltonian functions based on our tetrads have as many 
integrals of motion as degrees of freedom. Thus, it is 
shown that in the Schwarzschild limit these Hamiltonians 
describe an integrable system. We view this as an impor¬ 
tant test, as in general, for different tetrad fields the de¬ 
scription provides Hamiltonians non-equivalent be¬ 
yond the given approximation. In any Hamiltonian sys¬ 
tem the integrals of motion play a crucial role, when the 
integrability issue is studied. If we have several possible 
descriptions of the same system in the given approxima¬ 
tion, those respecting all background symmetries are the 
obvious choice. We use the case of spinning particle in 
Schwarzschild spacetime as such an exact problem with 
many integrals of motion to demonstrate shortcomings 
of certain coordinate-tetrad choices. Even though the 
considered approximations assume small spins, to clearly 
demonstrate (non-)integrability we also use large spin 
values in numerical tests. 

As for the Kerr spacetime, it was shown in [s^ that if 
the MP equations supplemented by the T SSC are lin¬ 
earized in the spin, an integral of motion associated with 
a Killing-Yano tensor appears. This led to the impres¬ 
sion that, up to linear order in the spin, the motion of a 
spinning particle is integrable in general [ 3 ^ . However, 
according to our numerical calculations, this seems not 
to be the case for the Hamiltonian function depending 
on the tetrad field choice introduced in [s^- 

This paper is organized as follows. In Sec. [2] we give a 


short overview of the Hamiltonian formalism introduced 
in [ 2 ^. After that, in Sec. [31 we present two different 
choices of coordinate systems, the Boyer-Lindquist and 
the Cartesian isotropic coordinates, for a tetrad corre¬ 
sponding to a ZAMO observer which is already given 
in [ 2 ^, [s^]. We analyze the properties of both with 
the help of analytical calculations and numerical integra¬ 
tions. Then we present our new tetrads in Kerr-Schild 
coordinates in Sec. [d] Finally, in Sec. |5l we summarize 
our results. We add a description of our numerical tools 
in the Appendix. 

We use geometric units, i.e., (G = c = I), and the 
signature of the metric is -l-,-|-). Greek letters denote 
the indices corresponding to spacetime (running from 0 
to 3), while Latin letters denote indices corresponding 
only to space (running from I to 3). 

2. THE HAMILTONIAN FORMALISM 

The Hamiltonian formalism in (2^ has been achieved 
by linearizing the MP equations of motion for the NW 
SSC. The MP equations describe the motion of a particle 
with mass satisfying the mass shell con¬ 

straint, and spin in a give n spacetime background 
Their reformulation in (40| reads 

£> 1 

^ = , ( 1 ) 

da 2 

D 

, ( 2 ) 

da 

where is the four-momentum, = dx^/da is the 
tangent vector to the worldline along which the particle 
moves, a is an evolution parameter along this worldline, 
and R^vkX is the Riemann tensor. The NW SSC reads 

= 0 , (3) 

where is a sum of timelike vectors. This sum in 
has the form 

Wy = p^ - m , (4) 

where e,y ^ is a timelike future oriented vector (through¬ 
out the article we use T instead of 0), which together 
with three spacelike vectors /, denoted by capital latin 
indices, is part of a tetrad field a- 

This tetrad field has to satisfy two conditions: the first 
condition ensures the orthonormality of the tetrad given 
by 

A P/ii/= PrA , (5) 

where prA is the metric of the flat spacetime and 
its analogon for the curved spacetime background. The 
capital indices are raised and lowered by the fiat metric 
PrA, the small ones by g^y. The second condition is 
implied by ([5]) and reads 

e>^Aey^=5^ 


( 6 ) 
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where 5^ is the Kronecker delta. 

When a tensor is projected on the tetrad field, then it is 
denoted with capital indices. For example, wa = 
is the projection of the time-like vector (|4]) on the tetrad 
field, i.e., 

wt = Pv g'" T -m , 

wj = Py g’' j . (7) 

Then, the spin tensor S^‘' projection reads 

g/ G,J . (8) 

In [i^, the authors do not work with this tensor but 
rather employ the spin three vector 

Si = 2'^iJL S^^ , (9) 


provides the contribution of the particle’s spin to the mo¬ 
tion, with 

+ E^ij^ , ( 18 ) 

and 

(Da = ojv S' a , 

(Dy = Pv - m it, , 

Pi = Pi , 

Po = -/ 3 ‘ p^ - aijrin? + , 

(Dt = Pu S' T -m , 

Qj = Pt, S' j . ( 19 ) 

The equations of motion for the canonical variables as 
a function of coordinate time t read 


where e/jL is the Levi-Civita symbol. 

Now, the Hamiltonian function H for a spinning par¬ 
ticle 


H = Hns + Hs , 
splits in two parts. The first. 


( 10 ) 


dx^ _ dH 
It 


dt 


dH 

dx^ 


dH 


dSi _ 
dt 


c 


( 20 ) 

( 21 ) 

( 22 ) 


Hns = d^P + a ^m2 + 7bp,p^ , (11) 


is the Hamiltonian for a nonspinning particle, where 


a 


/3* 
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gOO > 


gOrgOj 

^00 


( 12 ) 

(13) 

(14) 


and Pi are the canonical momenta conjugate to x® of 
the Hamiltonian COD. They can be calculated from the 
momenta pi with the help of the relation 


The phase space of a canonical Hamiltonian system 
is equipped with a binary operation, i.e., the Poisson 
bracket. If the dynamical system is subject to (sec¬ 
ondary) constraints ^i, the Poisson bracket has to be re¬ 
placed by the Dirac bracket [4lj | 

{Q,R}j,B-={Q,R}-{Q,d^}[c-^],^{dj,R} , (23) 

where Q and R are functions on phase space and C~^ is 
the inverse of the matrix consisting of the Poisson brack¬ 
ets of the set of constraints C = {^i, In the case of a 
spinning particle the constraints are given by the supple¬ 
mentary condition, here the NW SSC (d, and, in order 
to retain the symplectic structure, by the choice of the 
timelike body-fixed tetrad vector to be aligned with the 
four momentum 


P, = p. + E^taS^^ , 

= p, + (2E,tj— + Euc] Sl , (15) 

\ UJT ) 

where the spin connection 

Evta = — 2 (^9k\ S^ r + g'^ r Fk^a , (16) 

is a tensor which is antisymmetric in the last two indices, 
i.e., E„rA = —E^af, and Fki^a are the Christoffel sym¬ 
bols. The second part of the Hamiltonian, 


Hs 



a^vp^Ff \ 
sJrrS F Y^PiPj ) 


Sc , ( 17 ) 


Pm n 

X/. := e^T-= 0 , 

m 

where a is related to the local frame a by a Lorentz 
transformation (for more details see [l^). In order to de¬ 
rive the canonical structure of the phase space variables, 
the new defined momenta in eq. m are treated as 
functions of the kinematical momenta p^, the position, 
and of the spin which result in the following bracket re¬ 
lations 


{S^,S^}^^ = e^^^S^ + 0{S'^) . (24) 

All the other bracket relations between the variables van¬ 
ish at linear order in spin [^ . At this approximation. 
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even if the mass m? = —p^p^ is not a constant of motion 
for the exact MP equations with NW SSC, actually it 
scales quadratically in the particle’s spin (see, e.g., [1]), 
the mass is preserved at first order in the spin and treated 
as a constant in the linearized Hamiltonian formalism 



When we restrict the scheme to the linearized Hamil¬ 
tonian formalism, and consider the Pi no longer as func¬ 
tions but as independent phase space variables, then the 
terms of O (5^) are dropped in all the above Dirac brack¬ 
ets, i.e., in (EH) and all the other bracket relations be¬ 
tween the variables Profoundly, in the lin¬ 

earized Hamiltonian formalism a quantity / is a constant 
of motion, if it holds for its Dirac bracket with the Hamil¬ 
tonian function iJ 

{/,i?}^s=0 . ( 25 ) 

This means that if the system is evolved by the eqs. EOD- 
( 1221 ) . then the quantity I is preserved during the evolu¬ 
tion. 

The formulation provided up to this point is general, 
namely it does not depend on the specific coordinate sys¬ 
tem or on the specific tetrad field. These two factors, 
however, are essential for the Hamiltonian function dm. 
In particular, the nonspinning part of the Hamiltonian 
function m depends on the coordinate system which 
the metric is written in, while the spinning part (I17|) de¬ 
pends on the tetrad we choose. In Sec. [3] and Sec. |H 
we present three different combinations tetrad O coordi¬ 
nates for the Kerr spacetime background and discuss the 
advantages and shortcomings of the respective setups. 

3. THE HAMILTONIAN FUNCTION IN 
BOYER-LINDQUIST COORDINATES 

COMPARED WITH CARTESIAN ISOTROPIC 
COORDINATES 

3.1. A tetrad in Boyer-Lindquist coordinates 

A Hamiltonian function for the Kerr spacetime back- 

g und in Boyer-Lindquist (BL) has been provided in 
. The line element of the Kerr spacetime in BL coor¬ 
dinates reads 

ds^ = gtt df + 2 dt d(j) + dcf^ 

+ grr dr^ + gee dO^ , ( 26 ) 

with 

2Mr 

gtt — — 1 H— — , 

2aMr sin^ 6 

9tij> = r:; j 



S 

9rr — j 

gee = ^ , 


and 

S = cos^ 9 , 

A = — 2Mr , 

2 2,2 
w = r + a , 

A = — a^A sin^ 6* . (28) 

M denotes the mass and a the spin parameter of the 
central Kerr black hole. 

The tetrad field given in (2^ reads 



where for the small indices the numbers have been re¬ 
placed with the corresponding coordinates, i.e., t, r, 9, (j) 
stand for 0, 1, 2, 3, respectively. The proposed tetrad 
corresponds to an observer in the zero angular momen¬ 
tum frame (ZAMO) which intuitively yields a reasonable 
choice. Moreover, the coordinate system is based on the 
spherical coordinates in flat spacetime which respects the 
symmetries of the spacetime. In the Schwarzschild limit 
the above tetrad held reduces to (a —>■ 0) 

> 

, 

^ = r , 

= r sin 9 , (30) 

where /(r) = I — In the flat spacetime limit (M —>■ 
0 , a —>• 0) we get 



= r sin 9 . (31) 

This yields the flat spacetime in spherical coordinates. 

Let’s have a closer look at the dynamics in 
Schwarzschild spacetime. The corresponding metric in 
Schwarzschild spacetime results in 

ds^ = —/ (r) dt^ + f {r)~^ dr^ + {d9^ -f sin^ {9) d(j>^) , 

with / (r) = 1 — 2Mjr and the corresponding tetrad held 
is (1301) . The Hamiltonian 

H = Hns + Hs , 

is expressed in terms of the new phase space variables 
{r,9,(j), Pr, Pe, P^, Sf^) where stands for the spin 









5 


projected onto the spatial background tetrad in spheri¬ 
cal coordinates (reduced from the Boyer-Lindquist coor¬ 
dinates). All told, we have 


^ x/7M sin %) 

fix) ( C0S(6>) oBL _ J^oBL 

^ sin (0)^2 





(32) 




where Q = -p / (r) -f _^^p 2 _ 

In [ 2 ^ a criterion for the behavior of the Hamiltonian 
in the flat spacetime limit was introduced in order to 
check whether the choice of coordinates is a “good” one. 
Ideally, the contributions from the spin to the Hamilto¬ 
nian Hg vanish, since we no longer have curvature which 
the spin could couple to and the trajectory of the spin¬ 
ning particle should simply be the one of a straight line. 
Thus, the motion of the particle should be completely 
independent of the spin. However, in the case of spheri¬ 
cal coordinates the Hamiltonian is given by (15^ and the 
contribution from the spin part Hs does not vanish rep¬ 
resenting an evolution of the spin in the absence of spin- 
orbit coupling, as was noted in (^ . which might imply a 
coordinate effect for this choice of tetrad. Following the 
latter line of thought, we might say that the basis vectors 
are coordinate dependent, since they are oriented along 
the direction of the coordinate basis vectors in spheri¬ 
cal coordinates. Therefore, they introduce an additional 
evolution to the dynamical system which affects the equa¬ 
tions of motion for the spinning particle, i.e., the equa¬ 
tions of motion do not only contain the physical dynamics 
of the spinning object, but also the coordinate dynam¬ 
ics. On the other hand, the coordinate effect might not 
be the only interpretation, for instance for a long time 
the helical motion of a spinning particle with Pirani SSC 
in the flat spacetime was considered unnatural, until it 
was explained in terms of a hidden momentum in [ 13 . 
Anyhow, such effects make it harder to gain insights into 
the physical behavior of the particle’s motion, since it is 
not so easy to distinguish between coordinate effects and 
physical effects in the results. Therefore, we prefer to fo¬ 
cus on a more solid criterion for the Hamiltonian to check 
whether the choice of coordinates is a “good” one, and 
this criterion comes from the symmetries of the system. 

Generally, according to Noether’s theorem each space- 
time symmetry is related to a conserved quantity. In the 
case of spinning particles moving in a particular space- 
time geometry equipped with a symmetry described by 
a Killing vector the associated quantity conserved by 
MP equations reads 

C = ■ (33) 

In Schwarzschild spacetime we have three spatial Killing 
vectors yielding the three components of the total angular 


Figure 1: The left panel shows the relative error of AJ^, and 
the right of AJy as a function of time in logarithmic scale 
for the Schwarzschild background. The gray lines show the 
relative error of these quantities when the system is evolved 
using the Hamiltonian function corresponding to the tetrad 
dsni), while the black lines show the relative error of these 
quantities when the system is evolved using the respective MP 
equations. Both evolutions share the same initial conditions 
conditions, where a = 0, M = m = 1, and 5 = 1. 

momentum 

Jx=-ps sin {(j)) - cot (0) cos (</>) 

-1- gjjj ('5/^2 ^Qg 

— sin (^) , 

Jy =p 0 cos [(f)) — Pfj) cot (0) sin {(j)) -h sin (0)^ sin {(j)) 

+ rS^'' sin (0) cos (0) sin {<j)) + rS'"^ cos (</>) , 

Jz =P 4 > — f sin (0)^ (5'^’' — rS^'^ cot (0)) , 

where pi are the kinematical momenta and 5*-^ the spin 
components written in coordinate basis. In order to check 
whether the components of the total angular momentum 
are constants of motion within the Hamiltonian formu¬ 
lation we have to transform the expression to canonical 
variables Pi and Sf^ with the relations given in ([5]) and 
m- Therewith we obtain 

Jx = cos{(j)){Si^ csc(0) — P^ cot(0)) — Pe sin(0) , 

Jy = Pe cos(0) -I- sin((/))(5f ^ csc(0) — P^ cot(0)) , 

Jz = P<p , 

for the components of the total angular momentum, with 
which we may now compute the evolution equations for 
Ji via the Dirac brackets with the Hamiltonian. Indeed, 
they result in 

{jy,H)j,^=0{S^) , 

{Jz^ H}DB ~ 0 I 

and 
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Although we consistently keep the linearization in the 
Hamiltonian and the corresponding bracket structure, we 
find that the Dirac brackets for J^, Jy and Jz contain 
contributions from higher orders in the particle’s spin. 
Indeed, Jx and Jy start oscillating when the Hamilto¬ 
nian system corresponding to the tetrad field (l30ll is nu¬ 
merically evolved through the the equations of motion 
(OIIl) - (l^^ . It is visible from the relative error 

= i = x,y , (34) 

at time t of the Jx and Jy (gray line) in Fig. [U 
that the Hamiltonian function resulting from the tetrad 
(|30)l apparently violates the symmetry properties of the 
Schwarzschild spacetime. Consequently, the total angu¬ 
lar momentum J^ is not preserved, because the x and y 
components of the total angular momentum exhibit in¬ 
appropriate behavior. On the other hand, the respective 
evolution using the MP equation supplemented with NW 
SSC, instead, shows the expected preservation of the an¬ 
gular momentum components (black curves in Fig. [1}. 
This shows that even in the above linear in spin Hamil¬ 
tonian approximation a quantity is a constant of motion 
only when its Dirac brackets with the Hamiltonian are 
exactly zero, while when the brackets have contributions 
from the higher in spin orders, the quantities show no 
constancy. The violation of the expected symmetries re¬ 
sults in a system that exhibits chaotic motion (scattered 
dots in the left panel of Fig. HD, which contradicts with 
the integrability of the Hamiltonian for the spinning par¬ 
ticle on the Schwarzschild background we prove in sec¬ 
tion o It is true, however, that the relative error of 
the Jx, Jy components, and therefore of J^ scale with 
S (right panel of Fig. [2D. However, this should be antic¬ 
ipated since as S' —0 the system basically ignores the 
spin contribution and tends to reproduce geodesic trajec¬ 
tories. 

In this work we focus on the properties of the equa¬ 
tions of motion, and not so much on the astrophysical 
implication of these equations. The spin of the parti¬ 
cle makes the trajectories to deviate from their geodesic 
paths. Thus, we can interpret the spin as a perturba¬ 
tion parameter of the system. A constant of motion can¬ 
not dependent on the magnitude of the spin, even if the 
given value might be astrophysically irrelevant. This in¬ 
dependence from the spin magnitude holds also for the 
integrability of a spinning particle Hamiltonian (exclud¬ 
ing of course the case when S = 0). In our numerical 
calculations we measure the spin in units of m M, i.e. 
a = S/{m M), and set both masses to 1, thus the spin 
parameter cr = S' is dimensionless. Large values of the 
dimensionless spin, like S = I, might be astrophysically 
questionable, but do not have implications on the dy¬ 
namics (see, e.g., [il[i3 for relevant discussions). In 
our paper large values of the spin serve mainly as a tool 
to amplify the effects we want to point out, since these 
effects become less prominent when S' <C 1. 


or 



Figure 2: In the left panel is a detail from the surface of 
section 0 = tt/ 2, Pe > 0. The parameters of the orbits are 
H = 0.95, Jz = P,), = 3, S = M = m = 1, a = 0, the 
common initial conditions are cf) = 0, Pr = 0, Si = S 2 = 0, 
S = S 3 = 10~^, while the initial Py is defined numerically by 
the Hamiltonian function H. The right panel shows how the 
maximal value of the relativity error AJ^ = |1 “ d^(t)/J^(0)| 
for evolution intervals t = 10® scales with the spin measure 
S. 


Moreover, we find that 

{Ll + Ll + LlH)^^ = 0{S^) , 

with 

Lx = -pe sin{cj)) - cot {9) cos{cj)) , 

Ly = pg cos{cj)) - p^ cot{9) sin((/i) , 

Lz — P<j> , 

which also have to be rewritten in terms of the canonical 
momenta Pi. The conservation of the measure of the or¬ 
bital momentum of the linearized in spin MP equation 
in the case of the Schwarzschild spacetime background 
has been thoroughly discussed in for the Pirani SSC 
0. When the measure of the spin = SiS^ and the 
total angular momentum J^ are preserved, the integral 
of motion is equivalent to the conservation of L • S'. 
However, as in the case of the total angular momentum, 
we recover the same numerical problems for the measure 
of the orbital angular momentum. These two oscillations 
issues can be traced back to the coordinate dependence 
of the basis vectors in the spherical coordinate system, 
as we will see in the next subsection. 

So far, these coordinate effects have been investigated 
in Schwarzschild spacetime. Since the Schwarzschild 
spacetime is the nonrotating limit of the Kerr spacetime 
we would like to ensure that such coordinate effects can 
be eliminated in the nonrotating limit, i.e., the coordi¬ 
nate effects should vanish for nonrotating or slowly rotat¬ 
ing black holes. Thus, we were wondering whether there 
are more suitable choices of a coordinate system and of 
a tetrad for rotating black holes which do not show any 
unphysical coordinate effects in the Schwarzschild limit. 
Hence, the question arises as to which coordinates are 
best used? 
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Therefore, in the rest of Sec. [3] we study the Hamilto¬ 
nian formulation in an isotropic coordinate systems for 
the same kind of observer (ZAMO), introduced by [s^- 


The above relation between r and R holds outside the 
black hole’s horizon (5ll |. 

When we go back to the Schwarzschild spacetime 
where a —>■ 0, 


3.2. The Hamiltonian function in isotropic 
Cartesian coordinates 


A revised Hamiltonian function for the Kerr spacetime 
background in BL has been provided in (3^ . The for¬ 
mulation starts in Cartesian quasi-isotropic coordinates. 
The line element in these coordinates for an axisymmet- 
ric stationary metric is 


ds'^ = gtt dt^ 

+ 2 gtxdX dt + 2 gtydY dt + 2 gxydX dY 
+ gxx dX'^ + gyy dY'^ + gzz dZ"^ , (35) 

with 


9 tt 

9 tx 

9ty 

9xy 

9 xx 

9Yy 

9 zz 


e-^'' [B^cj^iX'^ + Y^) - 
e-^^uB^Y , 
-e-^''ujB^X , 

- e'^^^)XY 
A2Tr2 ’ 

A2 -hy2 ’ 

g2/.y2 ^2 

x^Ty^ ’ 

e2M 


(36) 


where w, e^, e'', B are functions of X, Y, Z. 

For this coordinate system the authors propose the 
tetrad field 


6/3 


T 


e/3 


1 


B ujY ^ e^^X^ + e-''B Y^ ^ 

{ei^ - e-'B)XY Y 

X2 + y^ /3 ’ 


ef) 


2 


+ 


BwX^, , {e>^-e-'B)XY 
ef^Y^ + e-^B X^ y 

-5 


e/3 3 = ef^sl , (37) 

corresponding to an infalling observer with zero 3- 
momentum. This tetrad becomes Cartesian, i.e., e^^ = 
1 , 6/3 in the flat spacetime limit. 

The Cartesian quasi-isotropic coordinates relate with 
the BL coordinate system through the transformation 


X = R (r) sin 6 cos (j) , 

Y = R{r) sin 6 sin (j) , 

Z = R{r) cos 6 , 

R{r) = ^{r-M + VA) . (38) 


ds^ =-fiR)dt^+ h{R){dX^+ dY^+ dZ^) , (39) 

the tetrad (l37l) reduces to the isotropic tetrad given in 

M- 


where 


e^T ^ = ./T(r)SI 


= y^^=(l + T17) > 


e/3' - + <^1 ’ 


M' 

mV 


(40) 


f{R) = 


r = i? 

(2i7-l) 


M 

XR 

2 


{2R+1) 
( XI 


2 ’ 


In order to check the behavior of these so called isotropic 
Cartesian coordinates (A, Y, Z) we analyze the conserva¬ 
tion of the constants of motion given by the symmetries 
of the system. The spherical symmetry of the spacetime 
can be described in Cartesian-like coordinates by the 
following three Killing vectors 

Using (I33p we thus get the three conserved components 
of the total angular momentum as a combination of the 
kinematical momentum and the components of the 
spin tensor On the other hand, in the canonical de¬ 
scription, the conservation of the total angular momen¬ 
tum 


JK = e^^^x^PM+SK , (42) 

is demonstrated by vanishing Dirac brackets 

{Jk,H}j,^ = 0 . (43) 

Contrary to the previous case the canonical momenta 
Pm and tetrad components of the spin appear in this 
formula. The relations between the two sets of quantities, 
the kinematical and the canonical, are given by m and 
([H]). By computing the difference of projection of (1551) 
and ([55)1 it can be shown, that if the Lie derivatives of 














the three spatial tetrad vectors obey the Cartesian-like 
rule 


L)e^iM = —^KLM and = 0 , (44) 

the two conserved quantities, one in kinematical variables 
and the other in canonical ones, are identical. Indeed, 
this formula holds in the flat Minkowski spacetime for 
the Cartesian tetrad = 8^, which naturally leads to 
the intuition, that a tetrad, that reduces to a Cartesian 
one in the flat spacetime, is a good tetrad choice. (In (l44l) 
the fact that the time component of the Killing vectors is 
required to vanish is explicitly stated, since it is written 
as a covariant, coordinate independent formula, but it 
has been derived using the above coordinate assumption.) 

The general condition can now be applied to the 
particular case of the Schwarzschild limit (1551) . As Lie 
derivatives can be written using partial rather than co¬ 
variant derivatives, one can easily check, that the tetrad 
field (HOI) satisfies (|541) . 

Yet, as an example, that the equivalence between total 
angular momentum expressed in kinematical and canoni¬ 
cal variables is not so obvious, let us consider a symmetry 
of the Schwarzschild spacetime w.r.t. rotation along the 
z-axis 


e = [0,-Y,Y,0] . (45) 

It yields the related component of the total angular mo¬ 
mentum 

J, =Xpy - Yp, + S^y (^h {R) + 

- z{xsy^ + Ys^^) . 

2R 

Here, pi represents the kinematical MP momenta and 
the coordinate spin components, the prime denotes 
a derivative with respect to R. In the Hamiltonian ap¬ 
proach we use the canonical momenta Pi and the pro¬ 
jected spin components , so it is necessary to perform 
a transformation from {pi.S"^^) to using the re¬ 

lations given in (fT51) and By doing this, terms pro¬ 
portional to h'{R) get absorbed into P^ and Py and the 
corresponding component of the total angular momen¬ 
tum can be written as 


J, = XPy-YP, + S3 . (46) 


The corresponding Hamiltonian in these coordinates, 
cf. reads 


H = Hns + Hs , 


with 

Hns 


Hs 


VfW) 


Vq 


( 1 + + m 


{l-s) 


(47) 

(48) 


and Q = m? + Notice, that setting M —>• 0, i.e., 

no gravitational field, we indeed obtain that the spin part 
of the Hamiltonian Hs becomes zero, as it should in the 
Minkowskian spacetime. 

Next, we can easily compute the evolution equations 
for the Ja;, Jy and Jz as 


{Jxi H} = 0 , 

H} gg = 0 , 

{Jz, H} gg = 0 , 

which is thus also true for the measure of the total angu¬ 
lar momentum . Moreover, it holds that = 0 

where = L^. + Ly + Li is the measure of the or¬ 
bital angular momentum. Its respective components 
are defined as Li = £ijkq^P^, with g* = {X,Y,Z) and 
H = {P,,Py,Pz). 

In fact, since the linearized in spin Hamiltonian 
system given by (071), 051 has five degrees of free¬ 
dom, the five independent and in involution constants 
of motion {Jz, , L^, , H) of the Schwarzschild limit 
make the system integrable. The integrability for the 
Schwarzschild background seems to result from the lin¬ 
earized in spin Hamiltonian approximation, because 
in 0 it has been shown that for the full MP equations 
with T SSC chaos appears for a spinning particle moving 
in the Schwarzschild background. However, the integra¬ 
bility seems to vanish in the Hamiltonian approximation 
once we turn on the spin of the central body. Namely, in 
the case of a Kerr spacetime background chaos appears 
again (scattered dots in Fig.[3|), which suggests the non- 
integrability of the corresponding Hamiltonian system. 
The existence of chaos in the Kerr background case for 
the Hamiltonian approximation is not just a confirmation 
of previous studies concerning the full MP equations with 
T SSC, see, e.g., [130. It shows that the linearized in 
spin Hamiltonian function given in is non-integrable. 
This result contradicts statements in the literature saying 
that up to the linear order in spin the motion of a spin¬ 
ning particle corresponds to an integrable system, see, 
e.g., [^. A thorough investigation of chaos in the Kerr 
spacetime for the linearized in spin Hamiltonian function 
will be provided in [4^. 

The above results match exactly the expectations we 
had from the symmetries. 

Up to this point, we have discussed the properties of 
a ZAMO tetrad in spherical and Cartesian coordinates 
in Schwarzschild spacetime. Taking the conservation of 
the constants of motion for numerical calculations as an 
important criterion to be satisfied, promising indicators 
for a “good” tetrad choice are the reduction to Carte¬ 
sian tetrad in flat spacetime as well as the vanishing of 
the spin dependent Hamiltonian. Two questions arise 
with this statement: First, are there other coordinates 
we may choose providing us with “good” tetrads, and 
second, since we were focusing on a ZAMO tetrad, we 
ask whether a non-ZAMO tetrad yields the same proper¬ 
ties if the coordinate basis is not changed. We expect the 
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and 


/ 


r 


2 M 


+ a/p^ + 4a^ 


-2 I -2 I -2 2 

X + 2 / + 2 — a 


(52) 


(53) 


For simplicity in the rest of the section we drop the bar 
notation over the KS coordinates. 

Independently on the tetrad field the choice of coordi¬ 
nates implies the non-spinning part of the Hamiltonian 

Hns = hPi + ay/rrp + PiPi - fa^{liPiY , (54) 


where hPi = S'^HiPj and 


r 

Figure 3: A detail from the surface of section 6 = n/2, Pe > 0. 
The parameters of the orbits are H = 0.9449111825230683, 
Jz = 3.5, S = M = m = l,a — 0.1, the common initial condi¬ 
tions are (j> = 0, Pr = 0, Si = 0, while by solving numerically 
the system Pe = —S2, Jz = P^ + S3, and S = + 5| we 

define the rest. 


1 

~7t+7 


(55) 


In the following we present two tetrad choices corre¬ 
sponding to different types of observers. 


4.1. ZAMO Tetrad 


properties of the tetrad to depend on the choice of the 
coordinates, so that in the following we take Kerr-Schild 
coordinates and analyze two tetrads, one ZAMO and one 
non-ZAMO tetrad. 


4. THE HAMILTONIAN FUNCTION IN 
KERR-SCHILD COORDINATES 


In the previous section, we focused on a tetrad field as¬ 
sociated to the observers with vanishing momenta Pi = 
0, i.e., zero angular momentum observers (ZAMO), in 
two different coordinate systems, isotropic Cartesian and 
Boyer-Lindquist coordinates. Therefore, it is reasonable 
to first consider such an observer in KS coordinates as 
well. Here, we choose a tetrad corresponding to an ob¬ 
server infalling with the radial velocity C/’’ = —af: 


The Kerr-Schild coordinates have the great advantage 
that they are horizon penetrating so that they are well 
behaved in the vicinity of the horizon, which simplifies 
numerical calculations in this domain, probably improv¬ 
ing the numerical treatment compared to isotropic co¬ 
ordinates for events in the strong field. Here we shall 
introduce a Hamiltonian function using the Kerr-Schild 
(KS) coordinates. The line element in KS coordinates 
reads [ij 


ds'^ = giivdx^dx’' , 
giif — ipti' f 1 1 / , 

where (0,1, 2, 3) correspond to {t, x, y, z), 

1 , 

r x + a y 

-b ’ 
f y — a X 

^2 _|_ q2 

z 
r 


k = - 


h = - 


(49) 


(50) 

(51) 


, 

e/3 ^ 1 - af) li 5°^ P [a~^- l) Ijlp . 

Again, this tetrad becomes Cartesian, i.e., ep’^ = Sj, 
ep ^ = S^p. in the flat spacetime limit, which is a first 
indicator for being a good tetrad choice. The next step 
is to analyze the behavior in the Schwarzschild limit a —>■ 
0. Then, following the procedure introduced in [ 2 ^, we 
obtain the Hamiltonian with 

= a(m-^r-P^ , (56) 

L-S 

J.3 ’ 

(57) 

where 


m=^m^+P^-^(r-Py , 

UJT = —m — Wi . (58) 


^Schw _ 


M 

m 


2a^ 


r® + 3a^ r • P 


— a — 
ujt 
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Since the total Hamiltonian is merely a function of 
certain scalar combinations of ^r, P, (where r = 

{x,y,z)), namely H = H (^r\'^,\P\'^,r ■ P,L ■ with 

Li = EijkT^P^ , we can deduce that 

\l + S,h] =0 , (59) 

t J DB 

by using the canonical structure of the variables. More¬ 
over, we would like to stress here again, that the conser¬ 
vation of in Schwarzschild spacetime is equivalent to 
the conservation oi L - S so that it suffices to express the 
Hamiltonian in terms ot L ■ S in order to show ([59|). In 
fact, it reflects the integrability of the system at linear 
order in spin. 

However, we cannot simply infer that J = L S is 
valid in the new canonical coordinates. The conserved 
total angular momentum is already given by (1331) and by 
the Killing vectors (IdTI) of the Schwarzschild spacetime 
we get 


J, = U + S, , (60) 

where the tilde denotes the quantities to be written in 
terms of the kinematical momenta pi {Li = Sijkr^p^) and 
the index i in Si refers to the coordinate basis. This 
relation is valid in KS coordinates, independent of the 
tetrad choice. 

In order to relate the conserved quantities to the 
canonical momenta Pi and the tetrad components of 
the spin S'/, we have to perform a transformation from 
{pi,Si) to {Pi, Si) using the relations given in (ITSl) and 
(O. Therewith, we indeed find the components Ji to be 
given by (gSl), which yields vanishing Dirac brackets for 
each component of the total angular momentum accord¬ 
ing to the argument mentioned in Section r3.2l In order to 
support this statement we performed a numerical check 
shown in Fig. |4l 




Figure 4: The left panel shows the relative error of Jx, and 
the right of Jy as a function of time in logarithmic scale for 
the ZAMO tetrad in KS coordinates as evolved by the Hamil¬ 
tonian with a = Q, M = m = 1, and S = 1. 


not accumulate over the integration time but stay at the 
same level. These results are similar to the ones obtained 
in isotropic Cartesian coordinates, so that the quality of 
the outputs is comparable. Therefore, if one can choose 
between KS and isotropic Cartesian coordinates, there is 
no preferred choice between those two in Schwarzschild 
spacetime. However, if the dynamics of plunging orbits 
is considered in a Kerr spacetime background, it may be 
more sensible to change to KS coordinates since they are 
horizon penetrating and avoid numerical divergence close 
to the horizon (see Appendix [B|) . 

Second, we consider the contribution from the spin 
part of the Hamiltonian in flat spacetime. From dSZD 
we see that for M —>■ 0 the contributions from Hs van¬ 
ish as it should. Hence, also additional coordinate effects 
which arise in spherical coordinates are avoided, further 
supporting such a choice of tetrad and coordinates. 


4.2. Non-ZAMO tetrad 


To simplify the Hamiltonian in Kerr-Schild coordinates 
we change to another tetrad field, which is not required to 
be a ZAMO observer. In particular, we take advantage of 
the fact that for certain observers no square roots appear 
due to normalization of the tetrad vectors 


en = 


e ^ - 




f f f f 

1 — — —/ —/ —/ 


/ 




^ 

f f f f ' 

— Li L] 1 1 _|_ £/ 7 L] 1 

2 27’ 2 « a;; ^ 

f f f f ' 

-^Iz, ^Izlx, ^Izly. l + ^lzlz 


(61) 

(62) 

(63) 

(64) 


where we use the definitions from above, cf. Eqs. dsni)- 
(15^ . This is the tetrad of an infalling ‘non-ZAMO’ ob¬ 
server, as the observer’s specific angular momentum 


e^T 


f±Y~ Ifa x^+y^ 

\d4>) 2 r r'^ + 




(65) 


and the observer’s radial coordinate velocity 


e T = 




Thus, we again compute the Hamiltonian in canonical 
coordinates up to linear order in spin given by 

H = P{j,fg -p Hso + Hss , ( 66 ) 

where Hus is given by (IMl) . 


It is immediately obvious that the conservation of these 
components is ensured up to numerical errors which do 


Hso 


af 


Mm — 2m{M — fr) 
2M m LOT 


r e^^^kpjS^ 


(67) 





















and 




11 


Hss = -- 


af 


( 68 ) 


dwr Mm (d^ll + r^) 1 

4 //z m ((m — affh)r + afhM) — 2M1^ (jnm + am^) 


S% 

'31 


+ 2 a [(M + 2 r)mflz + (3M — 2/r)P2] m 
+ 2 a (to + m) [MIz S^Pi — to (2 /r — 3M) 5^] 

- 2 (3 Mr - a^flz"^ - 3 fr^) \a{S^Py - S^P^) 


— {am + TO — afm) {S^ly — S^lx) 


Here, instead of (1551) . we used 


TO = v^to2 + PiPi - fa^{liPiY , 

TO = am — a^PiP , 

, f ~ frn\ 

lot = —m -1—TO , (69) 

a 2 


which together with the usage of components of P in¬ 
stead of coordinates significantly shortened expressions 
for Hso and Hss- All vector components are grouped 
in such a way that the relation {Lz + Sz,H}jjg = 0 is 
obvious. 

Again, the complete angular momentum conservation 
is restored in the Schwarzschild limit. Since Hns only 
depends on the chosen coordinate basis, it is still given 
by (1551) . The spinning part 


Hs = 


M ( M + 2r r-P 
TO I r (r -I- 2M) lot 


M 1-f 
wt1 + ^ 


L-S 

J.3 ’ 

(70) 


where to and lot are given by dMD, can again be writ¬ 
ten as a function H = H |Pp, r • P,L ■ so that 
we can follow the reasoning of the preceding subsection 
to obtain vanishing Dirac brackets (1551) . Therefore, we 
only have to check the equations for the components of 
the total angular momentum Ji in canonical coordinates 
{Pi, Si). Using the expressions for the total angular mo¬ 
mentum with respect to the coordinate basis dsni) , we 
again perform a transformation to the tetrad basis and 
the canonical momenta and recover relation (I46|) . Thus, 
in the Schwarzschild limit, the non-ZAMO tetrad in KS 
coordinates has the same numerical properties as the 
ZAMO tetrad, as expected, which is also visible in Fig. [5] 

Consequently, it seems to be a good choice of coordi¬ 
nate system for numerical investigations. 

It is of course also possible to rewrite the coefficients 
of the tetrad basis vectors in terms of any coordinate 
system without changing the general properties of the 
Hamiltonian system as long as the tetrad basis vectors 
remain oriented along the isotropic coordinate (Cartesian 
like) basis vectors. In (2^, it was already mentioned that 




Figure 5; The left panel shows the relative error of J^, and 
the right of Jy as a function of time in logarithmic scale for 
the non-ZAMO tetrad in KS coordinates as evolved by the 
Hamiltonian with a = 0, M = m = 1, and 5 = 1. 


the coordinate effects can be avoided by choosing the 
directions of the tetrad basis vectors along a Cartesian 
coordinate system. However, if the tetrad corresponds 
to a Cartesian frame, the spin variables remain Cartesian 
whereas the position and momentum variables are spher¬ 
ical ones. This approach is used in effective-one-body 
theory or post-Newtonian methods in order to compare 
the dynamical contributions from different orders in spin, 
(see e.g., [H, |4^) and may in fact also be used for the 
computation of the equations of motion from the Hamil¬ 
tonian. Nevertheless, in that case it is more sensible to be 
consistent in the choice of coordinates and spin variables 
so that the Dirac brackets can be used for the calculation 
of the equations of motion. This coordinate system does 
not necessarily adapt to the symmetries of the spacetime 
as we have seen. Generally, it is very useful to choose a 
coordinate system and corresponding basis vectors that 
do not imply coordinate effects if one aims at the analysis 
of the equations of motion. 


5. CONCLUSIONS 

In this work we have studied the Hamiltonian formal¬ 
ism of a spinning particle provided in with regard to 
numerical investigations of the equations of motion. It 
was already discussed in [2^ that this Hamiltonian for¬ 
malism does not only depend on the tetrad field one uses, 
but also on the coordinate system one chooses in order 
to express the tetrad field or the Hamiltonian function. 
Using the Dirac brackets to check the integrals of motion, 
we have shown that an unfortunate choice of the coordi¬ 
nate system can lead to a nonpreservation of quantities 
in numerical integration which should, according to the 
symmetries of the system and the linearized MP equa¬ 
tions, be conserved. However, we find that the type of 
the tetrads, i.e., whether the observer is ZAMO or fol¬ 
lows some other worldline which does not correspond to 
a ZAMO does not affect the general dynamical proper¬ 
ties of the constants of motion. In fact, we have exam¬ 
ined both kinds of tetrads and found no difference in 
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their ability to be numerically applied, i.e., they possess 
the same properties with respect to numerical computa¬ 
tions. However, the formulae for the spinning part of the 
Hamiltonian can be simplified and compactified, which 
we think is worthwhile to be mentioned. 

In order to obtain Hamiltonian systems without co¬ 
ordinate effects smearing the actual physical behavior 
in numerical solutions and still reliable in the vicinity 
of the central object’s horizon, two new horizon pene¬ 
trating Hamiltonian functions were introduced. Both of 
them were constructed on tetrad fields which were ex¬ 
pressed in Kerr-Schild coordinates. In particular, the 
non-ZAMO tetrad allows us to express the Hamiltonian 
both in Schwarzschild and Kerr spacetime in a simple 
and compact form. Future (numerical) work may profit 
from this explicit Hamiltonian. 

While studying the Dirac brackets in the Schwarzschild 
limit, we have shown that in this limit the Hamiltonian 
functions with acceptable properties are integrable. In 
particular, we have shown that the system’s five degrees 
of freedom admit five independent and in involution in¬ 
tegrals of motion. On the other hand, we have shown 
by a numerical example that chaos appears when we go 
to the Kerr background in the case of the Hamiltonian 
function proposed in [s^- This suggests that this lin¬ 
earized Hamiltonian function corresponds to a noninte- 
grable system in the case of the Kerr background. Indeed, 
the Hamiltonian formulation offers a wide range of appli¬ 
cations in the context of chaos and perturbation theory, 
such as Poincare sections or recurrence plots. In order 
to answer the question for chaos thoroughly, a detailed 
analysis of the motion of spinning particles in the Kerr 
spacetime described by the linearized in spin Hamilto¬ 
nian approximation is in progress 
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Appendix A: Numerical integration of the 
Hamiltonian equations 

Our numerical integrators rely on the considerations of 
our previous work, cf. [^, Secs. A and B], which can be 
extended to all Hamiltonians in this work. Let us briefly 
summarize the main points. 

All Hamiltonian equations of this work pos¬ 
sess a so-called Poisson structure, i.e., for 
y = {Pi,P 2 ,P 3 ,x^,x'^,x^,Si,S 2 ,S 3 )'^ € R®, they 


can be written as 


y = B{y)VH{y) , (Al) 


where B : R® —>■ R®^® is the skew-symmetric matrix¬ 
valued function 

/ 0 -/3x3 0 \ 

Biy) = /3X3 0 0 , (A2) 

V 0 0 B,{y)J 

with 

/I 0 0\ 

/3X3 = 0 1 0 , (A3) 

Vo 0 ij 

( 0 -53 52 \ 

Bi{y) = 53 0 -5i . (A4) 

V--52 5i 0 y 


D ue to this sp ecial structure, the spin length 5 = 
\/5i -b 5| -I- 5| is conserved along solutions of the equa¬ 
tions of motion ini. Thus, the three dimensional spin 
S = (5i, 52, 53 )^ can be represented by two variables a 
and ^ via 


S = 5 


/ a / 1 - cos(q) 
( \/l sin(a) 


(A5) 


see, e.g., One can then show, cf. [l^, that 




dH 

da 

dH 



(A6) 

(A7) 


holds. Hence, in the transformed variables z = 
(Pi, P2, P3, a:®, X®, a), the equations of motion take 

the symplectic form 


z = 
J = 


J-^VH{z) , 

/ 0 ^4x4^ 

^”-^ 4 x 4 0 J 


(A8) 

(A9) 


As a nice consequence, we can evolve the system with 
Gauss Runge-Kutta schemes which have already been 
shown to yield very good results for little computational 
costs in previous studies, see, e.g., fi^. 

An s-stage Gauss Runge-Kutta scheme is a collocation 
method, i.e., an implicit Runge-Kutta scheme 


yn+i = y™ + h^6i/(Y/ , (AlO) 

s 

Yi = y„-b/i^ay/(Yj) , i = l,...,s , (All) 
i=i 
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with coefficients 


II 

(A12) 

II 

(A13) 

where the stages ci,..., are chosen as 



(A14) 

with Ci being the roots of the Legendre-polynomial of de¬ 
gree s. Here, h denotes the time step size, Yi, i = 1,..., s, 
are the so-called inner stage values and y„ denotes the 
numerical approximation to the solution y at time r = 
nh. The functions li{t) are the Lagrange-polynomials of 
degree s. 

'.w=n^ ■ 

(A15) 


Gauss Runge-Kutta methods have a convergence order 
0 {h^^) which is the highest possible order among collo¬ 
cation schemes, e.g., [49| . Detailed information on their 
implementation is given in [13, Sec. 7], and [10, Chap¬ 
ters VIII.5 and VIII.6]. 

Very importantly. Gauss Runge-Kutta schemes almost 
exactly preserve the Hamiltonian throughout the numer¬ 
ical evolution, cf.[0. Fig. 10]. Furthermore, it is known 
from numerical analysis that the solution y^-i-i coincides 
with the value at t = h of the interpolation polynomial 
U(t) through the points (0,yn) and (ci, Yi)...(cs, Ys). 
This interpolation polynomial can be shown to stay 
0 {h^) close to the exact solution of the equations of mo¬ 
tion, see, e.g., [13, Chapter II. 1]. Therefore, we can 

conveniently calculate 0 {h'^) approximations to surface 
sections, such as the one presented in Fig. |3| above (for 
the procedure details see, e.g., na)- 

Appendix B: Plunging orbits 

The analytical properties of equations of motion have 
also an impact on the behavior of their numerical solu¬ 
tion. In Fig. ini we plot the relative error of the Hamil¬ 
tonian AH = |1 — H{t)/H{0)\ for plunging orbits using 
the Gauss Runge-Kutta method (Appendix]^ with fixed 
time step for equations of motion given by the Hamilto¬ 
nian function discussed in Sec. 13.21 (BL) (top plot) and 
the one discussed in section 14.21 (KS ) (bottom plot) as 


functions of the BL radius. The figure clearly shows that 
the KS case covers smoothly the horizon of the Kerr black 
hole (vertical dashed gray line), while the BL fails to do 
so by definition. 

The mass parameters are M = m = 1 and the Kerr 
parameter is a = 0.9. As initial conditions for the BL 
orbit we use r = 25, 0 = 7r/2, = 0, and ^ ^ ^ = 

0. The first two conditions for the velocities determine 


- 0.5 0 . 0.5 1 . 



Figure 6: The relative error of the Hamiltonian AH as a 
function of the BL radius r for two plunging orbits. The 
top plot shows the AH of a plunging orbit evolved by the 
Hamiltonian function discussed in section(BLl. while the 
bottom plot shows the AH of a plunging orbit evolved by the 
Hamiltonian function discussed in section IQ IKS'). In both 
cases we used masses M = m = 1 and Kerr parameter a = 
0.9. For more about the initial conditions refer to the text. 
The horizon is depicted in both plots by a vertical dashed 
gray line. 


the initial values of S 3 , and Pr- The remaining initial 
conditions are P^ = Sx = Sy = Pg = 0. The time step 
is AtsL = 0.1. The initial conditions for the KS orbit 
are x = 25,y = z = 0, and Px = Py = Pz = 0 with spin 
Sx = Sz = 0, Sz = 0.9. The time step is Atxs = 0.1. 

In fact, to test how much the pole-dipole approxima¬ 
tion of rotating body fails near the horizon has to be in¬ 
vestigated by other studies. But, one can assume that not 
having equations of motion singular at the horizon helps 
to obtain the right numerical results within the range of 
spins allowed by the linear-in-spin approximation. 
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